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ABSTRACT 

An online entry guidance algorithm has been 
developed using a predictor-corrector approach. The 
algorithm is designed for the Reusable Launch 
Vehicle (RLV) and is demonstrated by using the X- 
33 model. The objective of the design is to handle 
widely dispersed entry conditions and deliver the 
vehicle at the Terminal Area Energy Management 
(TAEM) interface box within an acceptable tolerance 
and without violating any of the vehicle physical 
constraints. Combination of several control variables 
is used in testing the performance and computational 
requirement of the algorithm. The control variables 
are the bank angle, angle-of-attack and the time for 
roll reversal. The bank angle and angle-of-attack 
profiles are the nominal profiles plus the 
perturbations in each direction. The initial guess of 
the bank profile is a 45 degrees bank angle with 
reversal at 360 seconds from liftoff. A six-element 
state vector is propagated to the TAEM interface box 
through the integration of the equations of motion 
(EOM). Altitude, heading and range errors are 
computed between the desired and the achieved state 
at the TAEM interface. These errors are used to 
correct the initial guess of the control variables. This 
process is repeated until the errors meet an acceptable 
level at the TAEM interface. Several numerical 
optimization methods are used to evaluate the 
convergent property of the predictor-predictor 
methodology. Successful results are demonstrated 
using the X-33 model. 


1. BACKGROUND 

An entry guidance algorithm based on the Space 
Shuttle drag-energy modulation approach has been 
developed for a RLV and demonstrated using the X- 
33 model. This guidance algorithm works well in 
terms of handling dispersions and successfully flying 
to TAEM interface while meeting the thermal, 
structural, and control power constraints. However, 
the thermal constraints for a metallic TPS (being used 
on X-33/RLV vehicle to reduce weight and improve 
operability) are different than those used for the 
Space Shuttle TPS. Also, the Shuttle-based guidance 
law uses gain-scheduled parameters, which requires 
off-line tuning for every mission. Since the gains are 
fixed for a flight, the guidance law only handles 
preplanned contingencies and lacks the adaptability to 
accommodate unforeseen circumstances during flight. 
The costs associated with this labor intensive 
approach can be substantial. 

The limitations of the existing guidance algorithm 
caused NASA MSFC to launch the Advanced 
Guidance and Control (AG&C) project [1] in April 
1999 with completion planned in March 2001. 
Different NASA centers and universities are 
developing several competing guidance algorithms. 
The goal is to determine, via simulation fly-off, which 
method has the highest potential on reducing the 
operational cost and increasing the mission safety of 
the RLV. Lockheed Martin Aeronautics - Palmdale 
(the Skunk Works) worked closely with NASA 
MSFC in defining the requirements for these 
algorithms. The candidate algorithms for ascent and 
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entry being considered are the hybrid ascent/abort 
guidance using optimal control and collocation 
techniques, ascent/abort guidance scheme based on 
Neighboring Optimal Control (NOC) techniques, 
trajectory following ascent guidance scheme based on 
continuous-time predictive control concept, predictor 
corrector entry guidance techniques, entry guidance 
by trajectory regulations about a reference trajectory 
and reduced order predictive/tracking entry guidance. 
This paper concerns the predictor-corrector entry 
guidance only. 

2. INTRODUCTION 

The goal of the second generation RLV program is to 
reduce the cost of placing the payload in orbit by an 
order of magnitude and the probability of loss of 
vehicle by two orders of magnitude relative to current 
Space Shuttle operations today. The objective of 
reducing the payload cost to $ 1000/lb and the loss of 
the vehicle to 1 in 10,000 is requiring the AG&C to 
be adaptive, robust, and reconfigurable. The 
proposed method in this paper will help to achieve 
the cost and reliability objectives. 

The predictor-corrector algorithm described in this 
paper is demonstrated by using the X-33 model. The 
X-33 is a subscale technology demonstrator for the 
RLV. A typical mission profile is shown in Figure 1. 
The predictor-corrector algorithm is applied to the 


descent portion down to the TAEM interface. The 
baseline entry guidance on the X-33 is a Shuttle based 
algorithm with modification to accommodate the 
peculiarity of the X-33 design especially in the 
metallic Thermal Protection System (TPS) area [2]. 

The major advantage of the predictor-corrector 
method is that the EOM are used explicitly in the 
formulation of the algorithm. Thus the vehicle can be 
modeled in as much detail as needed. Having the 
state vector available every integration step allows the 
calculation of the TPS thermal stress indicator and its 
use as an active constraint. Other areas of concern 
are the real-time implementation and the rate of 
convergence to feasible solution. 

Predictor-corrector algorithms have been evaluated 
for the past several years. Reference [3] describes 
numerical roll reversal predictor-corrector algorithm 
in aerocapture and precision landing guidance for the 
Mars Surveyor missions. Predictor-corrector 
algorithms have also been demonstrated for entry and 
landing missions. Reference [4] describes predictor- 
corrector scheme for a Space Station resupply/rescue 
vehicle. The remaining sections of this paper 
describe the general approach, algorithm 
development, optimization methods, and set of 
constraints, simulation results, conclusions, 
acknowledgements, and references. 
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Figure 1 - X-33 Mission Profile 
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3. APPROACH 

The attractive reason for designing entry guidance 
algorithm on line is the ability to include the thermal- 
stress indicator (TSI) as an active constraint. The TSI 
computation is a function of the temperature variation 
within the panels (AT) and the pressure variation (AP) 
across the TPS panel. This requires solving a 
simplified one-dimensional heat transfer equation 
through the panel. In addition* the TAEM interface 
box needs to be satisfied (i.e. altitude* heading, and 
down range.) Multiple steps approach is considered 
to start with simpler formulation to more elaborate 
scheme as we progress in achieving our goals. 

The first approach is concentrated on achieving the 
TAEM interface using the bank angle as the only 
control variable. A six-element state vector is 
propagated using transnational EOM to the TAEM 
interface defined by a desired box. The difference 
between the propagated state and the desired state at 
the TAEM interface is defined as the error state. 
Several iterations in modifying the bank angle profile 
are used to generate the functional relationship 
between the error and bank angle command. Several 
optimization schemes are tried to find the control 
variable solution using MATLAB simulation. 

The second step is introduced to allow for 
incorporating the TSI constraint. Additional control 
variables are needed to increase the degrees of 
freedom such that the TAEM interface and TSI 
constraints are met simultaneously. The nominal 
bank profile is described as a single doublet as shown 
in figure 2. The bank angle control variables are Agi, 
A g 2 and At changes for reversal time t r . The nominal 
angle-of-attach profile is described as a three level 
step function as shown in figure 3. The angle-of- 
attack control variables are Acq* Aa 2 , Aq, and At 2 . 


a 



a 



Figure 3 - a control variables 


4. ALGORITHMS DEVELOPMENT 

A key consideration to a successful predictor- 
corrector algorithm is establishing the sensitivity of 
the state error relative to the control variable and 
defining the performance index to solve for the 
control variable. Hundreds of trajectory solutions 
from entry to TAEM interface are generated to 
establish the functional relationships between the 
state error and control variables. Typical examples of 
heading error and down range error as function of 
Aai and A g 2 at a constant t r are shown in figures 4 
and 5, respectively. 



Figure 4 - Heading error as function of A G\ and Ag 2 
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Figure 5 - Range error as function of AO| and A Cz 

Several performance indices are considered in 
generating the control variables. They ranges from 
the straightforward approach of finding the zeros of 
the nonlinear function to minimizing a L2 norm of the 
combined errors. Consistent results are obtained 
using the later approach of minimizing the combined 
errors function/ with respect to the control variables, 
where: 
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A/? e is the down range error, 

SF { is the scale factor associated with AR C 
is the heading error 


SFi is the scale factor associated with AXc 
A/t e is the altitude error 

is the scale factor associated with Ah e 

Although it is recognized that real time 
implementation is dependent on the convergence rate 
of the solution, more emphasis is placed on the 
consistency of finding the solution. There are several 
ways that can speed the convergence of the solution 
by using specially tailored optimization schemes 
coded in C rather than MATLAB. 

5. SET OF CONSTRAINTS 

The TPS constraints during entry are the primary 
focus of guidance algorithms. The shuttle-based 
entry guidance [5] using a drag-energy modulation 
can accommodate stagnation temperature and 
reference heat rate as shown in figure 6. However, 
for metallic TPS panels the TSI constraint is the 
major concern. In order to implement that constraint, 
a one-dimensional heat transfer equation through the 
TPS panel is used. The input to that equation is the 
state vector evaluated at each time step of integrating 
the EOM. The output of that equation is the 
temperature differential across the TPS panel and the 
equilibrium wall temperature of the outer panel. Both 
these quantities are needed to calculate the TSI and 
the creep indicator (Cl) of the panel. The solution of 
the heat transfer equation is repeated at different 
panels that represent critical locations over the 
vehicle. Figure 7 show's the different body points 
used in calculating thermal constraints and the layout 
of the TPS metallic panels. The structural constraints 
are represented by the bending moment, shear, and 
torsion at the root of the canted fins and vertical fins. 
The load indicator is a measure of the severity of the 
structural loads are the normal acceleration (N z ) and 
the product of dynamic pressure multiplied by angle- 
of-attack (Qa). 

normal acceleration 



4 Figure 6- Entry Drag Limits 
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The 100-Series Body Points 


#112 - leeward centerline 



Figure 7- Aerothermal constraint body points and TPS layout 


The other set of significant constraints are defined at 
the TAEM interface box in terms of altitude, down 
range, and heading. 

The number of control variables describing the angle- 
of-attack profile and bank angle profile are more than 
or equal to the number of active constraints during the 
optimization. 

6. SIMULATION RESULTS 

A 3-DOF MATLAB simulation for the X-33 is used 
for the conceptual development of the predictor- 
corrector entry guidance. The simulation is validated 
against the MSFC test bed simulation called 
MAVERIC (Marshal Aerospace VEhicle 
Representation in C.) The propagation of the three 
transnational EOM is performed using the standard 
fourth order Runge-Kutta integration scheme. The 
optimization scheme to solve for the control variables 
subject to the set of constraints is the default option of 
MATLAB. Although it is recognized that the use of 
MATLAB simulation is not the most efficient way for 


real time computation, it was useful in determining 
the proof of concept. Several perturbed state cases 
for initial entry condition and different initial solution 
guesses are used to test the robustness of the 
predictor-corrector algorithm. 

The first set of results is obtained using bank 
modulation only to meet the TAEM interface 
constraints of downrange and heading requirements 
as shown in table 1. It was apparent that the altitude 
at TAEM interface was consistently on the low side. 
At that time the last step in the angle of attack profile 
(Aa 2 ) was introduced as an additional control 
variable. It w'as found that small change in Aa 2 (in 
the order of 2 degrees) was enough to bias the TAEM 
altitude to be inside the desirable box. The last set of 
results is obtained using both bank and angle-of- 
attack modulation to meet the TAEM interface and 
TSI constraints as shown in table 2. 

Figures 8 through 14 show' the time responses for the 
low energy case (V o -50O fps) using 4 control 
variables as summarized in Table 2. 
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Table 1 - Predictor-Corrector guidance with 3 control parameters, minimize L2 norm 



TAEM altitude, ft 


No. of iterations 
No. of function 
evaluations 


initial Acjj, deg 


initial Aa 2 , deg 


initial t,, sec 


final Aa t , deg 
final Aq : , deg 

final t r , sec 



Table 2 - Predictor-Corrector guidance with 4 control parameters, minize L2 norm 





initial Aa, deg 


initial sec 


final Act | , deg 


final Act 2 , deg 


final Aa, deg 


final tr, sec 
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Figure 8 - Angle-of-attack vs. time 
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Figure 9 - Bank angle vs. time 



Figure 10 - Altitude vs. time 
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Figure 1 1 - Relative velocity vs. time 
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Figure 12 - Range error vs. time 
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Figure 14 -TSI 105 vs. time 
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7. CONCLUSIONS 

The results obtained from MATLAB simulation 
demonstrated successfully the applicability of using 
the predictor-corrector entry algorithm for the X- 
33/RLV operations. The next step is to demonstrate 
how efficiently this algorithm can be coded and used 
in 6-DOF MAVERIC test bed including Monte Carlo 
simulation. Another key factor in adopting the 
algorithm to real time application is establishing the 
initial solution guess as function of the perturbed 
initial conditions. 
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